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Abstract 

The long-time behavior of the velocity autocorrelation function (VACF) is investigated by the 
molecular dynamics simulation of a two-dimensional system which has both a many-body interac- 
tion and a random potential. With strengthening the random potential by increasing the density of 
impurities, a crossover behavior of the VACF is observed from a positive tail, which is proportional 
to t^^, to a negative tail, proportional to — The latter tail exists even when the density of 
particles is the same order as the density of impurities. The behavior of the VACF in a nonequi- 
librium steady state is also studied. In the linear response regime the behavior is similar to that 
in the equilibrium state, whereas it changes drastically in the nonlinear response regime. 
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Transport coefficients are related to time integrals of the time correlation functions of ap- 
propriate quantities. Sufficiently fast decay of the correlation functions is therefore necessary 
for convergence of the transport coefficients, and it had been expected that the correlation 
functions are exponentially decaying functions in time. There are, however, certain systems 
in which correlation functions exhibit power-law decay at longtimeSjWhich is called the 



long-time tail 
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191]. In heat conduc- 



tion systems such as nonlinear lattices, hard-core particles, and Lennard- Jones systems, the 
long-time tails of the autocorrelation functions of the heat flux and the resultant system-size 
dependence of t.e t.e^a. eonduet.vitv have been observed nn,ne*any fl fl fl 1 

The long-time tails have also been observed for the velocity autocorrelation function 
(VACF) of a moving particle, which determines the self-diffusion coefficient in particle 



transport systems 



g particlej_v 

],y,y,y,y, 



10 



11 



12 



13 



15 



16 



17l . [isl . [l9j ] . One of such systems 



is the hard-core fluid where the long-time tail of the VAC 

t^^/^ (which we call the fluid-type tail) 0, y, y, y, y, \m 



is positive and proportional to 



ll| . Here d is the dimension of 



the system. Hence D is divergent in the two-dimensional fluid, which may be observed as 
logarithmic dependence of D on the system size. Another system is the Lorentz model, 
which describes the motion of a single point particle in a disordered system where immobile 
scatterers are randomly arrayed. The long-time tail of the VACF in the Lorentz model is 
3e negative and proportional to — (which we call the Lorentz-type tail) 
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171]. Although this tail does not lead to divergence of D even for d = 2, 



it might induce a system-size dependent term in D which decays only algebraically with 
increasing the system size. 

Although in the fluid system only a many-body interaction is present whereas in the 
Lorentz model only a random potential (which immobile scatterers generate) is present, there 
are many physical systems which have both of them. A typical one of such systems is the 
electric conduction system, where interacting particles correspond to electrons and a random 
potential is generated by impurities and defects. At room temperature, two-dimensional 
electron systems in semiconductors are well described as two-dimensional classical systems 
because the room temperature is smaller than the quantum level spacing in the confining 
(z) direction and larger than the Fermi energy in the x-y plane [2^ . It is widely known from 
experimental results that D of such systems is well-defined, independent of the sample size, 
to a good approximation. This fact suggests that the coexistence of a many-body interaction 
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(among electrons) and a random potential would change the long-time behavior of the VACF 
from the above-mentioned theoretical results which assume existence of only either one of 
them. 

In the present paper, using the model of a classical electric conduction system proposed 
by us {211, we compute the VACF in a two-dimensional system which has both a random 
potential (impurity scatterings) and a many-body [electron-electron (e-e)] interaction, and 
study how the long-time behavior of the VACF changes with varying the densities of impu- 
rities and electrons [22] . We also investigate the VACF in nonequilibrium steady states, and 
observe considerable changes in the nonlinear response regime. 

First, we analyze the equilibrium state. For this purpose we use a simplified version of 



our model 



2l[ |. It consists of two types of classical hard disks in a two-dimensional square 



box, the linear size of which is L. One type, which we call impurity, is immobile. Its radius 
and number density are denoted by Ri and n,, respectively. The other type, which we call 
electron, moves freely till it collides with another electron or with an impurity. The mass, 
radius, and number density of the electrons are denoted by nie, Re, and ng, respectively. 
We can change the strength of the random potential and the frequency of e-e collisions by 
varying the values of rii and rie, respectively. The model without the impurities corresponds 
to the hard-disk fluid whereas the model without the e-e interaction corresponds to the 
(non-overlapping) Lorentz model. 

Using the event-driven molecular dynamics (MD) simulation, we study the long-time 
behavior of the VACF, C{t) = {v(t) ■ v{0))/d, of an electron in an equilibrium state {d = 2). 
The data of the VACF shown below are computed as follows: First we calculate the VACF 
for each electrons and take the average over the electrons in the system for each configuration 
of the impurities. Then we average the results over five configurations of the impurities. The 
errorbars are the standard deviations among the results of the five configurations. 

We set rrie, Re, and the temperature k^T to be unity, and take L and Ri to be L = 480 and 
Ri = 0.4 in these units. The boundary conditions are set to be periodic in both directions. 
The impurities are distributed randomly under the restriction that the distances among 
them are larger than 2(i?e + Ri)- The initial positions of the electrons are so randomly 
arranged as not to overlap with the other disks, and their initial velocities are given by 
the Maxwell distribution of unit temperature. We calculate the VACF in the following two 
cases: (i) when we fix rig and change rij, i.e., we introduce the random potential into the 
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FIG. 1: Collision frequencies per unit time, K^-e (circles) and K^-i (squares), plotted against nj. 
The value of Tig is fixed to 0.0217. The sizes of most of the errorbars are smaller than the sizes of the 
symbols. The region in gray shows the crossover region between the fluid-type and Lorentz-type 
tails. 

system where only the e-e interaction is present, (ii) when we fix rii and change ng, i.e., we 
introduce the e-e interaction into the system where only the random potential is present. 

In case (i), we calculate the VACF of the system with n^. fixed to 0.0217 and rti varied 
from to 0.0434. In this case, the reduced electron density is ne{Re + ReY = 0.0868. In 
Fig. [1] we show the e-e and e-i collision frequencies (denoted by K^-e and K^-i, respectively). 
Here, Ke-e = (number of the e-e collisions per unit time) /rieL'^ and Ke-i = (number of 
the e-i collisions per unit time)/neL^. From this result the duration time Xg-e (= 1/Ke-e) 
of the e-e collisions is evaluated to be about 12 when = 0. The numerical results of 
the VACF are shown in Fig. [2l We observe that C(0) ^ 1 = k^T/rrie, implying the law 
of equipartition of energy, and that the VACF decays exponentially at initial times and 
then decays algebraically at longer times. There exists the positive long-time tail for rii < 
2.17 X 10~^. The exponent is about —1 at Ui = 0, which reproduces the result for the hard- 



disk fluid 



0, y, 0, y, y, 0, 



ll| . With increasing rii this fluid-type tail becomes weaker (i.e., 
the decay becomes faster). For rii ^ 4.34 x 10~^ it disappears, and instead, the negative 
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FIG. 2: The VACF plotted against t, for various values of nj. The parameters are the same as 
those in Fig. [TJ Top: A semilog plot around initial times. Inset of the top: A magnified plot on 
linear scales. Bottom: A double logarithmic plot of the absolute value of the VACF. In the bottom, 
for clarity, the data of neighboring parameters have been shifted along the time axis by a factor 
of two. The dashed and dotted lines are guides to the eye, which are proportional to t^^ and t^^, 
respectively. 
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FIG. 3: Collision frequencies per unit time, K^-e (circles) and Ke-i (squares), plotted against ng. 
The value of rii is fixed to 0.0434. The sizes of most of the errorbars are smaller than the sizes of 
the symbols. The crossover region would exist at larger values of Ue than the range shown here. 



long-time tail appears. The VAC 
the Lorentz model 1^, Q, 
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decays as —t ^ for > 10 ^, which is the same result as 



id llTj l . This observation suggests that the Lorentz- type 



tail exists even when T^-e is the same order as Te-i = 1/Ke-i, the duration time of the e-i 
collisions. 

To see this fact more clearly, we compute the VACF for case (ii), with Ui fixed to 0.0434 
and rig varied from 4.34 x 10~^ (i.e., a single electron in the system, corresponding to the 
system without the e-e interaction) to 0.0434. In this case, the reduced impurity density is 
ni{Re + RiY = 0.0851, and K^-e and K^-i are shown in Fig. [31 From this result r^-i ~ 5 when 
the e-e interaction is absent. We show the numerical results of the VACF in Fig. HI It is 
seen that at long times the VACF decays as — at ne = 4.34 x 10~^, which corresponds to 
the Lorentz model, and that this Lorentz- type tail also exists for larger values of n^. Thus 
the Lorentz-type tail appears even when t^. 

Il6l . |l7j, the origin of the 



^ Te-i. (see Fig. ^ 

According to the phenomenological approaches [g, |9|, llO, 111 



long-time tails is the coupling of hydrodynamic modes in the systems. In the fluid the 
conservation of both the momentum and particle number contributes to the hydrodynamic 
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FIG. 4: A double logarithmic plot of the absolute value of the VACF, plotted against t, for various 
values of rig. The parameters are the same as those in Fig. [3l The red line (ng = 4.34 x 10^^) 
represents the case of no e-e interaction, which corresponds to the Lorentz model. For clarity, 
the data of neighboring parameters have been shifted along the time axis by a factor of two. The 
dotted lines are guides to the eye, which are proportional to 



modes Id, \1 



In the Lorentz model, by contrast, only the particle-number conservation 



contributes to the mode, which couples with the static mode generated by the configuration 



of the impurities pj7|. Although the absence of the interaction between particles (the e-e 
interaction) was assumed in ref. 17, our results suggest that a similar mode-coupling theory 
would be valid whether or not the interaction between particles is present if a system has 
the dynamical mode from particle-number conservation and some static modes couple with 
it. 

According to this interpretation, the crossover of the long-time behavior from the fiuid- 
type tail to the Lorentz-type one with increasing rii (Fig. [2]) would be explained as follows. 
To treat the electron system as a fiuid, it is necessary that the e-e collisions occur more 
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frequently than the e-i colhsions. That is, Xe-j > McTe-e is a necessary condition for re- 
garding the velocity field to be hydro dynamic, where Mc (> 1) is a crossover ratio of the 
collision frequencies, which would be almost independent of the densities and the system 
size. Therefore, when Tg-j/re-e ^ M^. both types of the long-time tails are present, and we 
observe the fluid-type one manifestly because it is stronger than the Lorentz-type one. With 
increasing nj, Te-i becomes smaller, and when Te-i/Te-e ^ the fluid-type tail disappears 
because the hydrodynamic description of the velocity field is not valid. The Lorentz-type 
tail, by contrast, survives even in this regime due to the particle- number conservation. 

In Fig. [2] the crossover takes place around 2.17 x 10~^ ^ f^i ^ 8-68 x 10^'^. Thus from the 
corresponding values of Tg-j and r^-e in Fig. [H we evaluate Mg to satisfy 

2 < Mg < 10. (1) 

We can translate Mg into the crossover value of the density ratio as 

{ne/ni)c = Mc(Te-i/(Te-e (2) 

if we assume the approximate relations r^-i ~ l/veUiae-i and Xg-g ~ l/fgngCig-g. Here, fg is 
the thermal velocity of an electron, and (jg-j and (jg-g are the cross sections of the e-i and e-e 
scatterings, respectively. Since Oe-ijoe-e = {Re + Ri)/2Re for hard-disk interactions, eq. 
is translated as 2 < {ne/ni)c < 7. These values are consistent with Fig. [2] (where rzg is fixed 
to 2.17 X 10~^), which shows the vahdity of eq. ([2]). In Fig. H {rii = 4.34 x 10"^), Ue/rii 
ranges from 10^^ to 1 in our calculation, which is smaller than {rie/nijc. Therefore, the fact 
that we have observed only the Lorentz-type tail in Fig. H] is also consistent with the above 
interpretation. 

Finally, let us examine the long-time tail of the VACF in nonequilibrium steady states. 



In order to realize a nonequilibrium steady state, we use our original model [2l|]. That is, 
in addition to the electrons and impurities there exist another type of moving disks (which 
we call phonons) in a rectangular box of size x Ly. The mass and radius of the phonons 
are rup = 1.0 and Rp = 1.0, respectively. In the x-direction we apply an electric field E to 
the electrons to induce the electric current, and the boundary condition in this direction is 
periodic. The boundaries in the y-direction are potential walls for the electrons, and thermal 
walls for the phonons, which simulate the heat transfer outside the system (conductor) to 
keep it in a steady state. The interaction between any pair of the disks is taken to be 
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FIG. 5: Left graphs: The mean velocities of an electron versus electric fields. The broken lines 
represent the slopes of the linear response coefficients. Right graphs: The absolute values of the 
VACF for steady states under various values of the electric field, plotted against t. The data in 
red, green and blue are the results for the equilibrium state, the linear response regime and the 
nonlinear response regime, respectively. The dashed and dotted lines are guides to the eye, which 
are proportional to and respectively. The number densities of the impurities are (a) rii = 
[ne/rii > {ne/ni)c], (b) = 0.00213 [ue/rii ~ {ne/ni)c], and (c) Ui = 0.0267 [ne/rij < {ne/ni)c]. 
Unlike Figs. [2]and[31 we have not shifted any of the data along the horizontal axes. 
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Hertzian. The equations of motion are solved by the Gear's predictor-corrector method. As 
shown in ref. 21 and in the left graphs in Fig. [5] (in which the charge of an electron is set 
to be unity), this model has well-defined nonequilibrium steady states even in the nonlinear 
response regime. We calculate the VACF defined by 

^<'^^' = ((„.-K»^> 

in a nonequilibrium steady state, where is the component of electron's velocity parallel to 
E. We set = 750 and Ly = 125, and fix the number densities of the electrons and phonons 
to be He = 0.0160 and Up = 0.00533, respectively. From the results for the equilibrium state 
{E = 0), we find that the crossover of the VACF with increasing is observed even when 
the phonons are present. Then we compute and compare the VACFs in the equilibrium 
state, the linear response regime and the nonlinear response regime, for various values of 
rie/rii. The right graphs in Fig. [5] depict the results. It is seen that in the linear response 
regime the behavior of the VACF is almost the same as in the equilibrium state. In the 
nonlinear response regime, by contrast, we observe drastic changes in the behavior. For 
rie/rii > {ne/ni)c [Figs. [5](a) and [5](b)] the fluid-type tail is enhanced as the electric field 
becomes larger. For Ue/ni ~ {ne/ni)c [Fig. [5](b)] in particular, although we cannot observe 
the fluid-type tail clearly either in the equilibrium state or in the linear response regime, it 
appears manifestly in the nonlinear response regime. For ne/rii < {ne/ni)c [Fig. M^c)] the 
Lorentz-type tail appears from shorter times than in the equilibrium state. We think that 



this is partly because t^-i becomes shorter with increasing E \21\\ . 

In summary, using the MD simulation we have studied the long-time behavior of the 
VACF for interacting electrons in a two-dimensional disordered system. For the equilibrium 
state we have found that with increasing the density of the impurities the long-time tail 
crosses over from the positive tail proportional to t^^ (the fluid-type tail) to the negative 
tail proportional to — (the Lorentz-type tail). The Lorentz-type tail survives even when 
we increase the density of the electrons up to the same order as that of the impurities. These 
results imply that in our two-dimensional system the diffusion constant does not diverge in 



the limit of infinite system size [23|. There may remain, however, a small term, which 
depends weakly (algebraically) on the system size. We have also found that the crossover 
occurs when r^-i/Te-e ~ or Ue/rii ~ M^ae-i/ (ye-e, where 2 < < 10. 

For nonequilibrium steady states, although in the linear response regime the VACF be- 
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haves similarly to that in the equilibrium state, its behavior changes considerably in the 
nonlinear response regime. There the fluid-type tail is enhanced for smaller rii whereas the 
Lorentz-type tail appears from earlier times for larger rij. 

We expect that these summarized results would be insensitive to details of the systems, 
and that elaborate experiments and a kinetic theory with two-parameter and Ui) density 
expansion would support our results. Furthermore, our estimate of Mc would be useful for 
analyzing not only the long tails but also other physical phenomena of many-particle systems 
because one can judge whether the system can be described as a fluid. 
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